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Abstract. This paper treats two topics: matrices with sign patterns and Jacobians of 
certain mappings. The main topic is counting the number of plus and minus coefficients in 
the determinant expansion of sign patterns and of these Jacobians. The paper is motivated 
by an approach to chemical networks initiated by Craciun and Feinberg. We also give a 
graph-theoretic test for determining when the Jacobian of a chemical reaction dynamics has 
a sign pattern. 



1. Introduction 

This paper treats two topics: matrices with sign patterns and Jacobians of certain map- 
pings. The main topic is counting the number of plus and minus coefficients in their de- 
terminant expansion, but other types of results occur along the way. It is motivated by an 
approach to chemical networks initiated by Craciun and Feinberg, see [CF05,CF06], and 
extensions observed in [CHWprept]. 

1.1. Determinants of Sign Patterns. The first topic, see ^ is purely matrix theoretic 
and generalizes the classical theory of sign definite matrices [BS95]. This subject considers 
classes of matrices having a fixed sign pattern (two matrices are in a given class iff each of 
their entries has the same sign (or is 0)), then one studies determinants. Call a sign pattern 
a matrix A with entries which are ±^ij or 0, where Aij are free variables. To a matrix B 
we can associate its sign pattern A = SP{B) with ±Aij or in the correct locations. If 
A is square, then the determinant of A is a polynomial in variables A^j, which we call 
the determinant expansion of A. We call a square invertible matrix sign-nonsingular 
(SNS) if every term in the determinant expansion of its sign pattern has the same sign. 
There is a complete and satisfying theory of these which associates a digraph to a square 
sign pattern and a test which determines precisely if the matrix is SNS, see [BS95]. 

In this paper we analyze square sign patterns and give a graph-theoretic test to count the 
number of positive and negative signs in their determinant expansions; Theorem 12.91 We 
extend the result to nonsquare matrices and call our test on a matrix the det sign test. 
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1.2. Jacobians of reaction form differential equations. The second topic, ^ in this 
paper apphes this to systems of ordinary differential equations which act on the nonnegative 
orthant M^q in 

(1.1) l^/w, 

where / : M>g ^ M°'. The differential equations we address are of a special form found in 
chemical reaction kinetics: 

(1.2) I - S.(.), 

where S" is a real dx d' matrix and ?; is a column vector consisting of d' real- valued functions. 
We say that system (11. ip has reaction form provided it is represented as in (11.20 with 

v{x) = {vi, . . .,Vd') and 

(1.3) Vj depends exactly on variables Xj for which Sij < 0. 

Call S the stoichiometric matrix and the entries of v{x) the fluxes. We always assume 
the fluxes are continuously differentiable. 



Our second main result, Theorem 13.21 describes which S have the property that the 
Jacobian matrix /'(x) = Sv'{x) has a sign pattern, meaning that each entry flj{x) has sign 
independent of x in the orthant. The characterization is graph-theoretic, clean and elegant. 
The question was motivated by works of Sontag and collaborators [AnS03, ArS06, ArS07]. 

Our third main result here. Theorem 13.151 when specialized to square invertible S counts 
the number of plus and minus signs in the determinant expansion of the Jacobian f'{x) = 
Sv'(x) of a reaction form f{x) = Sv{x) in the terms of a bipartite graph associated to S and 
the det sign test. We use this to obtain results on the determinant expansion for general 
nonsquare 5*. 

We present many examples which illustrate features of our results and limitations on how 
far one can go beyond them. 



1.3. Chemistry. The reaction form differential equations subsume chemical reactions where 
no chemical appears on both sides of a reaction, e.g. catalysts. Furthermore, in many 
situations all fluxes Vj{x) are monotone nondecreasing in each Xi when the other variables 
are fixed, that is, v'{x) has all entries nonnegative. This happens in classical mass action 
kinetics or for Michaelis-Menten-Hill type fluxes. See [Pa06] for an exposition. 

A key issue with reaction form equations is how many equilibria do they have in the strict 
positive orthant ffi>o- observed in [CF05,CF06,CF06iee] that in many simple chemical 

reactions the determinant of /' has constant sign on the positive orthant and as a conse- 
quence of a strong version of this, any equilibrium which exists is unique. Other approaches 
exploiting this determinant hypothesis (under weaker assumptions) are in [BDB07, CHW- 
prept]. Roughly speaking, if the determinant of the Jacobian /' does not change sign on 
a compact region Q, then degree theory applies and bears effectively on this issue; the full 
orthant M.^^ can easily be approximated by expanding f2's. 

The degree argument is very flexible and probably extends to many situations. Fragile, 
however, are establishing constraints on the sign of the determinant. A key tool is the 
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determinant expansion of (11.21) . namely, the expression det(S'y(x)) as a polynomial in 
the functions Vij{x), which are the entries of the matrix function V{x) = v'{x). The main 
issue is the sign of the terms in the determinant expansion, are they all the same or if not are 
there few "anomalous" signs. In [CHWprept] it is observed that in each example of Craciun 
and Feinberg the determinant expansion has very few anomalous signs. When this happens, 
then it gave some methods one could use to prove existence and uniqueness of equilibria. 
For example, if the determinant expansion has one minus sign and many plus signs, and if 
the monotonicity condition Vij{x) > holds, then det(S'\^(a:)) is positive on large regions 
(which in particular situations can be estimated). 

Our main results. Theorem 13. 151 etc.. on det(S'y(a;)) were motivated by a desire to develop 
tools for counting anomalous signs. While the paper is not aimed at chemical applications, 
many of the examples of matrices 5* we use to illustrate our work are stoichiometric matrices 
for chemical reactions. 

The paper [BDB07] identified chemical reaction determinant expansions initiated by [CF05, 
CF06] with classical matrix determinant expansion theory and sign patterns. This is de- 
scribed in the book [BS95] and pursued into new directions in a variety of recent papers such 
as [BJS98,CJ06,KOSD07]. The bipartite graph conventions in this paper are a bit different 
than conventional, but were chosen to be reasonably consistent with [CF06]. 

The authors wish to thank Vitaly Katsnelson for diligent reading and suggestions. 

2. Matrices with sign patterns 

This section gives the set-up and our main results on sign patterns as described in the 
introduction. 

Let t{A), respectively 'm±(A), denote the number of terms, respectively ± signs, in the 
determinant expansion of the square sign pattern A. Recall a sign definite (SD) matrix 
A is one with either m-{A) = or m_(.(/l) = or det(/l) = 0. The number of anomalous 
signs m{A) of a square sign pattern A is defined to be 

mlA) := min{m_{A),m_^.{A)}. 

We say A is j-sign definite if it has j anomalous signs, that is m{A) = j. 

Lemma 2.1. m{A) = m{A) and m^{A) = m_(A) if A is obtained from A by: 
interchange of two rows, or interchange of two columns, or multiplying a row by a minus 
sign, or multiplying a column by a minus sign. 

Proof. Obvious. ■ 

Question (J-sign): Given a sign pattern S we are interested in whether every square 
submatrix is sign definite or more generally in getting an upper bound J on the j for which 
S contains a j-sign definite square matrix. 

We shall settle this question and give an even more refined result for square sign patterns 
A which counts m±{A). 
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2.1. Basics on graphs, matrices and determinants. To this end we revert to graphs. 
Given a sign pattern S let G{S) denote its signed bipartite graph. This is a simplified 
version of the species-reaction (SR) graph in [CF06]. It is a signed bipartite graph with one 
set of vertices C{S) based on columns and the other set of vertices R{S) based on rows. 
There is an edge joining column c and row r iff the (r, c) entry Src of S is nonzero. The sign 
of this edge is the sign of Src- If two edges meeting at the same column have the same sign, 
they are called a c-pair. By a cycle we mean a closed (simple) path, with no other repeated 
vertices than the starting and ending vertices (sometimes also called a simple cycle, circuit, 
circle, or polygon). A cycle that contains an even (respectively odd) number of c-pairs is 
called an e-cycle (respectively o-cycle). Recall a matching in a bipartite graph is a set 
of edges without common vertices. Equivalently it is an injective mapping from one of the 
vertex sets to the other. A matching is called perfect if it covers all vertices in the smaller 
of the two vertex sets. A k x k square submatrix A of S corresponds to k column nodes 
C{A) and k row nodes R{A); there is an associated sub-bipartite graph G{A) of G{S). 

Example 2.2. The following is an example taken from [CF05, Table 1.1. (v)] which illustrates 
these definitions. Given 



S 



-1 -1 



-1 




1 







-1 





1 





1 

-1 
-1 










-1 
1 

-1 






the signed bipartite graph G{S) is as follows: 



C3 



CI 



C2 



C4 



Here the dashed lines denote positive edges and full lines represent negative edges. 

The edges C3— R4 and C3— R3 are a c-pair, while C3— R4 and C3— R2 are not a c-pair. 
The cycle C3— R4— C4— R3— C3 has one c-pair, so is an o-cycle. On the other hand, the cycle 
CI— Rl— C2— R3— C3— R2— CI has two c-pairs and so is an e-cycle. ■ 



We use repeatedly the basic fact of linear algebra that if A is an n x n square matrix, then 

n 

(2.1) detA=Y, sign(f^) JJ ^^-(0 ' 



cr&S„ 



1=1 
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where Sn is the group of permutations on {1,2, ■ ■ ■ , n} and Aij denotes the term of 
A). 

Lemma 2.3. The bipartite graph G{A) of a square sign pattern A has no perfect matching 
iff det A = 0. 

Proof, det A ^ iS the determinant expansion will have at least one nonzero term, say 
yli^o-(i) ■ ■ ■ An,a{n) for some a G Sn- So all the v4j^o-(j) are nonzero, hence row 1 is connected to 
column cr(l), . . ., row n is connected to column cr(n) and this yields a perfect matching for 
G{A). m 

Remark 2.4. The same argument shows that: 

(1) If the bipartite graph of a rectangular matrix does not have a perfect matching, then 
the determinants of all of its maximal square submatrices are 0. 

(2) The number of terms in the determinant expansion of a square sign pattern A is the 
number of perfect matchings of G{A). ■ 

Remark 2.5. Note the following: 

(1) Without loss of generality for (J-sign) we can remove the second, third etc. or any 
colinear column from S. Thus there are no colinear columns in S. This is true 
because, (a) if A contains linearly dependent columns then det A is 0. (b) if A is 
the same as A except one column is removed and replaced (in any order) by a scalar 
multiple of that column then the only possible change in m±. is m±{A) = mzf{A). 

(2) Any cycle in G{S) can be embedded in a square submatrix of S. ■ 

2.2. SNS matrices vs. e-cycles. Let us call a square invertible matrix sign-nonsingular 
(SNS) if every term in the determinant expansion of its sign pattern has the same sign [BS95, 
Lemma 1.2.4]. If all square submatrices of (a not necessarily square matrix) S are either SNS 
or singular, then S is strongly sign-determined (SSD). 

Proposition 2.6. A sign pattern S is SSD iff the signed bipartite graph G{S) has no e-cycle. 

Proof. This fact is essentially classical, cf. [BS95, Theorem 3.2.1]. Also, it is a special case 
of Theorem 12. 9[ ■ 

2.3. Many Cycles: Square Matrices. Next we turn to the more general situation where 
e-cycles occur in G{A). The bipartite graph G{A) enables us to count the number of positive, 
negative and anomalous signs in the determinant expansion of a square sign pattern A. We 
permute and re-sign to make sure all diagonal entries of A are negative. Thus these diagonal 
entries correspond precisely to a perfect matching W in G{A). A cycle in G{A) that contains 
each edge (c, VV(c)) in G{A) corresponding to any column c it touches, is called interlacing 
with respect to W, or W-interlacing for short. 

Remark 2.7. To a given signed bipartite graph G we can associate (uniquely up to trans- 
position and a permutation of rows and columns) a sign pattern A with G{A) = G. The 
number of anomalous signs of a signed bipartite graph G with equipollent vertex sets is 
defined to be m{G) = m{A). 
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Example 2.8. Consider the following two graphs. 



(g) 



C3 



I 



C2 



@ 



@ 

I 

I 
I 



C3 



CI 



C2 



Graph Gi 



Graph G2 



CI 



Graph Gi admits only one perfect matching W. Namely the set of edges {CI— R2, C2— Rl, 
C3— R3}. Hence its only cycle Rl— C3— R2— C2— Rl is not W-interlacing. The sign pattern 
associated to Gi is 

Bu 
B = —B21 —B22 




■^13 

B23 
-B33 



As det{B) = -B12B21B33, m(G'i) = m{B) = 0. 



Graph G2 on the other hand admits three perfect matchings. For instance, with respect to 
the matching {C1-R3, C2-R2, C3-R1}, the cycle R1-C3-R2-C2-R1 is W-interlacing, 
while the cycle CI— R3— C3— R2— C2— Rl— CI is not. The sign pattern associated to G2 is 








Cl2 


-Ci3 


c = 


—C21 


—C22 


C23 












Since det(C) = -C12C23C31 - C12C21C33 + C13C22C33, m{G2) = m{C) = 1. 

Note that the number of interlacing cycles depends on the matching chosen. For instance, 
the graph G3 



C2 



C3 



C4 



CI 



R2 
Graph G3 

with the matching {CI— R3, C2— Rl, C3— R2, C4— R4} admits three interlacing cycles, 
while it has four cycles interlacing with respect to the matching {CI— R2, C2— Rl, C3— R3, 
C4-R4}. ■ 



The following theorem gives our det sign test counting the number of signs in the deter- 
minant expansion of a square sign pattern A in terms of G{A). For the sake of simplicity it 
is stated for matrices with nonzero diagonal entries. This causes no loss of generality since 
such a matrix can be obtained from any square invertible matrix with a permutation of rows. 

Theorem 2.9. Let A he a square sign pattern with nonzero diagonal elements. The diagonal 
gives us a perfect matching W that is fixed. 

(1) The number of terms, t{A), in the determinant expansion of A is one plus the cardi- 
nality of the set of all sets of disjoint yV -interlacing cycles of G{A). 
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(2) Let e be the sign of the product of the diagonal elements of A. Then the number of 
terms of sign — e in the determinant expansion of A, m_^{A), equals the cardinality 
of the set of all sets of disjoint W -interlacing cycles that contain an odd number of 
W -interlacing e-cycles. 

Remark 2.10. By disjoint cycles we mean cycles with no common vertices. The empty set is 
not counted set of cycles. ■ 

Remark 2.11. As observed in Example 12. 8[ the number of interlacing cycles depends on the 
matching W chosen. However, the numbers t{A), m±{A) and m{A) obtained from Theorem 
12.91 are (clearly) independent of W. ■ 

The special case of Theorem 12.91 where m{A) = is settled by [BS95, Theorem 3.2.1] 
which is due to Bassett, Maybee and Quirk [BMQ68]. 

The count of the signs in the determinant expansion is simple in extreme cases, as the 
following corollary shows. 

Corollary 2.12. Let A be a square sign pattern with nonzero diagonal entries. The diagonal 
induces a perfect matching W. Let e denote the sign of the product of the diagonal elements 



(1) Suppose that G{A) has t cycles interlacing with respect to W and each pair of cycles 
has a nonempty intersection. Then the number of terms in the determinant expansion 
of A is l-\-t and m-^{A) is the number of W -interlacing e-cycles. 

(2) Suppose that there are t > 1 cycles of G{A) each of which is W -interlacing and all 
are pairwise disjoint. Then the number of terms in the determinant expansion of 
A is 2* and the number of anomalous signs is either {if all W -interlacing cycles 
are o-cycles) or 2^~^. In the former case, = 2* and in the latter case m_{A) = 



Proof. For (1) note that every set of disjoint interlacing cycles contains only one cycle. (2) 
By Theorem 12.91 (1). the number of terms in the determinant expansion of A is just the 
number of all subsets of {1, . . . , t}, i.e., 2*. 

For the second part of the claim we will compute the number m-^{A). Let r be the 
number of e-cycles among the t interlacing cycles. Of course, if r = 0, then there will be no 
anomalous signs. So assume r > 0. There are t — r interlacing o-cycles. Since the cycles are 
pairwise disjoint, we have by Theorem 12.91 that a set consisting of some of the t interlacing 
cycles contributes a term with sign — e to the determinant expansion of A iff it contains an 
odd number of the r e-cycles. Thus to find m-e{A) we multiply the number of ways we can 
choose an odd number of e-cycles from the r e-cycles by the number of ways we can choose 
any number of o-cycles from the t — r o-cycles. The number of ways we can choose an odd 
number of e-cycles from the r e-cycles is 



of A. 



m+{A) = 2 
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To simplify this, notice that = (—1 + 1)^' = ^^=q (^,)(— l)'^ imphes 

L%^J / r \ 1 ^ 



ir— 1 



k=l ^ ' k=0 

The number of ways we can choose a subset of o-cycles from the t — r o-cycles is 2*~^. Thus, 
m_,(A) = 2^-1 ■ 2^-'' = 2*-i and hence m{A) = m±{A) = 2*"^ ■ 

Example 2.13. We now show how to determine when the determinant expansion of a square 
sign pattern A has no or one anomalous sign. Let us assume that all diagonal entries of A 
are nonzero and thus induce a perfect matching W. By Theorem 12.91 (2). m{A) = iff G{A) 
contains no W-interlacing e-cycles. 

We claim that m{A) = 1 iff G{A) contains exactly one W-interlacing e-cycle and no W- 
interlacing cycles disjoint from it. Clearly, (<^=) follows from Theorem 12.91 For the converse, 
note that if G{A) contains at least two W-interlacing e-cycles, then m{A) > 2 by Theorem 
12.91 (2). Similarly we exclude the possibility of only one W-interlacing e-cycle with other 
W-interlacing cycles disjoint from it. ■ 



2.3.1. Proof of Theorem l2.9i As preparation for the proof of the theorem, we briefly recall 
some well-known facts about Sn- A S'„-cycle s = (si ■ ■ ■ Sm) is the permutation mapping 

Si I— i> S2 I— ^ ■ ■ ■ I— ^ Sm, I— > Si 

and fixing {1, . . . , n} \{si, ■ • ■ Sm} pointwise. To avoid collision with cycles in various graphs 
appearing in the paper, we call these cycles S'„-cycles. 

Example 2.14. For instance, the S'4-cycle a = (12 4) is the mapping 

/I 2 3 4 \ 



V 2 



1 / 



while the mapping 



7 \ 



7 / 



can be written as (12 4) (3 5). 



Every permutation a & Sn can be written uniquely (up to the ordering in the product) as a 
product of disjoint S'„-cycles. Conversely, every set of disjoint S'n-cycles gives a permutation 
in Sn- 

Lemma 2.15. If a = ri ■ ■ ■ Tm is a factorization of a & Sn into disjoint Sn-cycles, and 
Ti = {TiiTi2 ■■■ TitJ, then 

n m ij 

sign(cT)j]A,.w=n(-i)*'"'n^-.--.-+i n ^^^^^ 

i=l j=l i=l k<^{nj} 

[with the convention Tji^+i = Tji)- 
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Proof. If a contributes to the determinant, then every Tj induces a cycle of G{A). For 
instance, the cycle G'(Tj) corresponding to Tj is defined to be the subgraph 

R{ni) - C{n2) - R{n2) - C{n^) i?(T,iJ - C(r,i) - R{m) 

of G{A). The other ingredient is sign(cr) = sign(ri) • ■ ■sign(rm) and sign(7:j) = (—1)*'"-'^. ■ 

Proof of Theorem \2.9[ Let A he n x n. Statement (1) follows from Remark 12.41 To see why 
(2) is is true, we invoke the determinant expansion formula (12. ip . For convenience we assume 
that all diagonal entries of A are negative. Then e = (—1)" and we count the number of terms 
with sign — e. Each term x = sign(cr) nr=i ^i,o-(«) expansion gives us a permutation 

cr e Sn- Since every permutation can be written uniquely as a product of disjoint S'„-cycles, 
we obtain a set of disjoint cycles ri, . . . , with a = ri ■ ■ ■ r^. Say Ti = {th Ti2 ■ ■ ■ tuJ. By 
the previous lemma, 

x=ri(-i)*-^n^.. ...... n 

Observe that the sign of a product of the form YIILi Ar^-^rji+i equals 

(^_l)nu™ticr of c-pairs in G(ri) _. g[gi^(^Q (^j-.'j'j 

Taking into account that all diagonal entries are negative, the sign of x then equals the sign 
of 

e 

(_l)ELi(i.-i)(_i)-ELii. JJsign(G'(r,)). 

i=l 

This simplifies further to 

e 

1=1 

In order for the term x to have sign (—1)"^^, (— 1)^ ]^^^-^ sign(G'(rj)) must not be equal 1. 
We will show this is the case iff the number of e-cycles among tq, . . . , is odd. 

Case (1): Suppose i is odd. Then 

e e 
{-lYl[sign{G{r,)) = -l ^ l[sign{G{n)) = 1 

1=1 1=1 

,^ ^ j^^#(o-cyclcs among to,...,T() 

<^=^ # (o-cycles among Tq, . . . , r^) is even 
<^=^ # (e-cycles among tq, . . . ,Te) is odd, 

since i is odd and # (e-cycles) = ^ — # (o-cycles). 
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Case (2): Suppose I is even. Then 

s, ^_-|^^#(o-cycles among TO,.. .,Tf) _ _-|^ 

# (o-cycles among tq, . . . , r^) is odd 

# (e-cycles among tq, . . . , r^) is odd, 

since £ is even and # (e-cycles) = £ — # (o-cycles). ■ 



2.4. Many Cycles: Nonsquare Matrices. The graph-theoretic test described in §2.31 
gives a det sign test settling Question (J-sign). In this section we extend the det sign test 
to nonsquare sign patterns S. 

A cycle has the property that the number of rows it touches is the same as the number 
of columns it touches. A set of cycles is called balanced if the number of rows they touch 
is the same as the number of columns they touch. Every balanced set of cycles picks out 
a square submatrix A oi S and hence induces a sub-bipartite graph G{A) of G'(S'). Such 
a submatrix and the sub-bipartite graph are both said to be balanced. Note each column 
and row of A appears in at least one cycle in G(v4). 

Proposition 2.16. For every square invertible submatrix B of a sign pattern S there is a 
balanced square submatrix A of S with m{A) = m{B). In fact, A can be chosen to be a 
submatrix of B. 



Proof. Suppose B is the smallest square submatrix of S violating the conclusion of the 
lemma. After permuting rows we assume B has nonzero entries on the diagonal. Since B is 
not balanced, either a row or a column of B does not appear in any cycle in G{B). Without 
loss of generality we assume this to be row 1. 

Since we assume that row 1 does not appear in any cycle of G{B), for a & Sn with 
a = Ti- ■ - Tg, where are disjoint S'„-cycles, the corresponding term in the determinant 
expansion x = sign(cr) nr=i -^«,o-(i) will be zero if 1 appears in one of the r^. Hence the 
nonzero terms x will correspond to permutations a with cr(l) = 1. In other words, Bi i will 
get picked from row one. So by removing row and column one from B we obtain a smaller 
matrix Bq with m(So) = rn{B). By the minimality assumption on B, there is a balanced 
square submatrix A of Bq with m{A) = m[Bo) = m{B), a contradiction. ■ 



By this proposition, the answer J to Question (J-sign) equals the maximal number of 
anomalous signs obtainable from a balanced square submatrix of the sign pattern S. So the 
algorithm for finding the desired upper bound J is as follows. Consider sets of balanced 
cycles in G{S). Each of these induces a square submatrix A of S. If G{A) admits no perfect 
matching, we continue with another set of balanced cycles. Otherwise we count the number 
of anomalous signs in det A by the procedure described in Theorem 12.91 of §2.3[ The highest 
possible count obtained is the desired sharp upper bound J. 
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3. Reaction form differential equations and the Jacobians 



Now we turn to studying systems of reaction form (RF) ordinary differential equations 
whicli act on the nonnegative orthant M>g in M'^: 

(1 T 

(3.1) ^ = /(x) = Sv{x), 

where / : M>q — S* is a real d x d' matrix and v is a. column vector consisting of d' 
real- valued functions. 

The differential equation (I3.ip has weak reaction form (wRF) provided V{x) := v'{x) 
satisfies Sij > ^ =0. If a differential equation has wRF, then it has reaction form 

provided Sij = ^ Vji{x) = and Sij < ^ Vji{x) ^ 0. The flux vector v{x) is monotone 
nondecreasing (respectively, monotone increasing) if ^{x) is either for all x E ]R>q 
or nonnegative (respectively, positive) for all x G M^'q. 

This section analyzes two properties the Jacobian of f{x) might have. First we say ex- 
actly when f'{x) has a sign pattern (Theorem 13.21 and Corollary 13. ip . Secondly we give a 
method based on §2] for counting the number of plus and minus coefficients in its determinant 
expansion. 



3.1. Sign pattern of the Jacobian. We first say precisely when the Jacobian of a reaction 
form dynamics respects a sign pattern and find that it does surprisingly often. 

Corollary 3.1. Given a reaction form differential equation 

Tt = ^"^^^ 

with monotone increasing flux vector v{x). The Jacobian Sv'{x) respects the same sign 
pattern for all x G M"q if the bipartite graph G{S) does not contain a cycle of length four 
with three negative edges. Conversely, if G{S) does contain such a cycle, then some matrix 
S arbitrarily close to S, possibly S itself, produces Sv'{x) which fails to respect the same sign 
pattern for all x in the orthant. 



The corollary is an immediate consequence of Theorem 13.21 which operates at a higher 
level of generality and requires the definition we now introduce. 

Here and in the sequel, U will denote the flux pattern assigned to S. It is a d' x d matrix 
with each entry being or a free variable Uij] the {i,j)th entry of f/ is iff Sji > 0. In case 
the differential equation (13.11) satisfies RF and the flux vector v{x) is monotone increasing, 
U is the sign pattern of V{x). 

Theorem 3.2. Let S be a real d x d' matrix and U the corresponding flux pattern. 

(1) The differential equation (13.11) has wRF iff each diagonal term in SU is a negative 
linear combination of monomials in Uij . 

(2) SU of a wRF differential equation admits a sign pattern {that is, each entry of SU is 
a positive or negative linear combination of monomials in Uij) whenever the matrix 
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(3.2) 



S does not contain a 2 x 2 submatrix with the same sign pattern as 



+ 1 
-1,0 



-1 
-1,0 



or 



-1 

-1,0 



+1 
-1,0 



or 



-1,0 -1,0 
+ 1 -1 



or 



-1,0 
-1 



-1,0 

+1 



" +1 


-1 " 




' -1 


+ 1 " 




' -1 


-1 " 




' -1 


-1 " 




-1 


or 


-1 




or 


_ +1 




or 


-1 


+ 1 _ 


-1 


-1 


-1 



Here —1,0 stands for either —1 or 0. 

(3) SU of a RF differential equation admits a sign pattern iff the matrix S does not 
contain a 2 x 2 submatrix with the same sign pattern as 

(3.3) 

Equivalently, in terms of the bipartite graph, G{S) does not contain a cycle of length 
four with three negative edges. 

(4) The entry {SU)ij is nonzero iff there is some k with Sik 7^ and Sjk < 0. If SU 
admits a sign pattern, then sign{{SU)ij) = sign(5'ifc). 

Proof (1) Write S = 5"+ — S'„ for real matrices S^, with nonnegative coefficients satisfying 
the comphmentarity property {S+)ij{S^)ij = 0. Diagonal entries of SU are of the form 
SijUji which meets the negative coefficient condition iff '^j{S+)ijUji = iff {S+)ijUji = 
for all i,j. This uses that the Uj are free variables, so no cancellation can occur. Thus 
{S)ij > iff {S+)ij 7^ implies Uji = which is the wRF condition. 

(2) The (2, j)th entry of SU does not have a sign pattern iff {S+U)ij 7^ and {S^U)ij 7^ 0. 
{S+U)ij = J2ki^+)ikUkj, so {S+U)ij 7^ iff for some k, {S)ik > and Ukj 7^ 0, i.e., {S)ik > 
and by wRF Sjk ^ 0. Similarly, {S^U)ij 7^ iff there is some i with {S^)ie 7^ and Uej 7^ 
so Sj£ ^ 0. Taken together this implies that the 2x2 submatrix of 5* given by rows i,j and 
columns k, i has the same sign pattern as one of the matrices in (13.21) . 

(3) This follows as in (2) by using that if RF holds, then Ukj 7^ iff Sjk < 0. Also Uij 7^ 
iff Sji < 0. 

(4) From {SU)ij = J2k ^ikUkj it follows that {SU)ij 7^ iff there is k with Sik 7^ and 
Ukj 7^ 0. Due to the construction of U, Ukj 7^ iff Sjk < 0. This proves the ffist part of the 
statement and the second follows immediately since Ukj is positive. ■ 

Our next step is to introduce several different types of determinant expansions. 

3.2. The core and other determinant expansions. For S e R'^'"^' with r := rank(5) 
we define the core determinant to be 



(3.4) 



cd(^) := lim det{SU -tl). 



Let i? be a matrix whose range is the orthogonal complement of the range of S. We also 
use the formula 

det{SU-BB^) 



(3.5) 



Cod{S) :-- 



det(i?i?^|(i^s)^)' 

which does not depend on which B we select and equals cd(S') (see Proposition 13. 3p . so we 
call both the core determinant. The Craciun-Feinberg determinant expansion [CF05] 
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is defined to be 

dd{S) := det{SU - tl) with t fixed, e.g. t = 1. 

For more details on the relationship between the Craciun-Feinberg determinant expansion 
cfd(S') and the core determinant cd(S') we refer the reader to §11 

Proposition 3.3. Let C G W^^'^' , let D be a d' x d matrix with possibly symbolic entries and 
suppose C has rank r. Define 

_ det{CD- BB^) 

"^"det(i?S^|(i^c)x)' 
where B is a matrix whose range is the orthogonal complement of the range of C . Then 

(1) a is independent of which matrix B whose range is the orthogonal complement of the 
range of C is used to define it. 

(2) a is the the determinant of the compression of CD to the range of C . 
(3) 

a = lim -J— detiCD - tl). 
Proof. We consider all matrices in the basis imC _L (imC)"*". Then 



CD 



C-D|imC * 





and BB^ 



* 
BB^ 



(imC)J 



CD - BB^ 



CD-tl 



Hence 

C-D|imC * 

-B-B^|(imC)-L 

This implies a = det(C-D|im5) and is thus independent of B. For (3), observe that 

{CD-tI)Uc 

— ^-^l(imC)^ 

As the size of the second diagonal block is {d — r) x [d — r), 

^ detiCD-tl) = det{{CD-tI)\imc) 



t^- 

Sending t — yields the desired conclusion. ■ 

Definition 3.4. A column c in iS G M'^^'^' is called reversible if — c is also a column of S. 
(Many matrices coming from chemical reactions have reversible columns.) We call — c the 
reverse of c. 

Example 3.5. Let us consider an example which is a slight modification of [CF05, Table 
l.l.(i)]: 



C3,C4 



C5 



ci,C2 — (m) 
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The corresponding stoichiometric matrix S and the vector v{x) are given by the following: 



an 


-On 








-Ol3 




' hxl'' 


021 


-021 


022 


-022 







k2x\^'xf^ 








O32 


-032 


033 


v{x) = 


hxT 


-041 


041 





















-052 


052 







hxl'' 



s 



Note some of the columns of S are reversible. This phenomenon is captured in the graph 
by listing two columns that are reverses of each other in a common rectangular box. For 
example, C3 and C4 appear in the same box and in fact columns 3 and 4 are reverses of each 
other. Sign of 5*33 is the same as the sign of 5*23 and both appear in the graph as a solid line. 
5*53 has sign opposite to these and so appears in the graph as a dashed line. This is also true 
for C4. Other dashed vs. solid lines of the graph coming from a box with reversible columns 
follow the same pattern. 



The corresponding V{x) and U are as follows: 



V{x) 



x 



,111- 



X 



ai3 





X2^^ Oil ^2 










X 



X2^ ^(12X^2 











1-4 



X 



«22— 1^132 



^022^4 



^122 ^132 



"^032^4 





























U21 


U22 























f/35 





f/42 


f/43 








f/51 



























^052 /cs 








u 



For generic choices of the numbers a^j the matrix 5" will be of rank 3 and this is what we 
focus on. A straightforward computation gives 

cd(S') = -2ai3a4ia52[/i4^735f^5i - 2013021052^22^/35^51 - 2013022041^714^42^/51 

-2ai3a3204i[/i4Kl3?751 - 2ai302i032[/22f7i3^751 + 2011022033^/22^/43^^51 

Hence there is potentially one anomalous sign in cd(iS). However, 

-2013021032^/22^/13^/51 + 2aiia22033[/22Ki3?/51 = 2(aii022 033 - Oi302l032)[/22f/43^51 

SO cd(5') has one, respectively no anomalous sign, depending on whether 011022033 — 013021032 
is positive, respectively nonpositive. ■ 



Example 3.6. Suppose S 



. Then SU — I admits a sign pattern; it is a 2 x 2 



—On On 
— O21 O21 

matrix with all entries negative. Hence the determinant expansion of its sign pattern has 
one anomalous sign by the det sign test. However, cfd(iS') and cd(S') have no anomalous 
signs. ■ 
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3.3. Formulas for determinants of products of matrices. For a matrix A, A{a\6) will 
refer to the submatrix of A with rows indexed by a and columns indexed by 6. 

Recall the Binet-Cauchy formula for the determinant of the product AB of a m x n 
matrix A and a. n x m matrix B: 

(3.6) det{AB) = J2 det(A(all|5)) det(B(5|all)). 

5C{l,...,n} 

I S I — m 

(If m > n, then there is no admissible set S and the determinant det{AB) is zero.) 

Combining Proposition 13.31 with the Binet-Cauchy formula we obtain 
Lemma 3.7. For S E W''^^' having rank r, the core determinant is given by 

(3.7) cd{S) = i-lf-'' det{S{a\(3)) det{U{(3\a)). 

\a\,\(3\=r 

Proof. Use (jM)) and [ PQ GH ] = [ P G] 

(3.8) det{SU - U) = ^ det ([ ^ -tl ] {al\\6)) det 

\S\=d 

where d is the number of rows of S. Since rank S is r, f^-'' factors out of det ( [ 5 -tl ] (all|5)) , 
so limt_^o det ([ S —tl ] (all|(5)) exists. Let us look at terms of degree d — r in t in (13. 8p . 
det ([ S —tl ] (all|(5)) will be of degree d — r in t ifS 6 will consists of exactly r columns (3 
of 5*. If a denotes the set of rows of S that do not hit any of the columns of —tl chosen by 
f3, then 

det{[S -tl ] (all|(5)) = i-lY-''t'^-'^det{S{a\P)). 

It is clear that such pairs (a, (3) are in a bijective correspondence with all S that pick r 
columns of S. Hence 

det{SU - tl) = i-ty-'' J2 det{S{a\(3)) det{U{/3\a)) + (higher order terms in t). 

\a\,\f3\=r 

Dividing by t'^~'^ and sending t ^ proves (13. 7p . ■ 

Formulas (13. 7p and (13. Sp are in contrast to cfd(S') which is given by the more complicated 
expression 

r 

(3.9) cfd(5) = 2]] {-ty~" det{S{a\f])) det{U{i3\a)) 

s=l \a\ = \l3\=s 

The fact is known (cf. [CF05], [BDB07, proof of Theorem 4.4]) and its proof follows the line 
of the proof of Lemma 13. 7[ 

For the chemical interpretation of the core determinant vs. the Craciun-Feinberg determi- 
nant see our ^ 

Lemma 3.8. The number of anomalous signs in cd(S') is at most the number of anomalous 
signs in cfd(S'). 



Q 

H 



to get 



U 
I 



(5|all) 
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Proof. By looking at the formulas (13.71) and (13.91) it is clear that each term appearing in 
cd(S') also appears (multiplied with f^"^) in cfd(>S'). Terms w from cd(S') have degree r in 
the Uij^s. All terms in cfd(iS') not coming from terms in cd(>5') have degree < r in the t/j^'s. 
Thus there is no cancellation and the statement follows. ■ 

Remark 3.9. Example 13. 211 shows that the number of anomalous signs in cd(S') can be strictly 
smaller than the number of anomalous signs in cfd(5'). 

If B is the sign pattern associated to the graph Gi of Example 12. 8[ and S = [ B —B ] , 
then cd(S') has no anomalous signs, whereas cfd(S') has one anomalous sign. We leave this 
as an exercise for the interested reader. ■ 



3.4. Generic matrices and the reduced iS-matrix. In this section we introduce some 
basic definitions and illustrate them with an example. 

Definition 3.10. A matrix A is called weakly generic if its rank r is maximal among all 
matrices with the same sign pattern. If, in addition, all r x r submatrices of A are weakly 
generic, then A is called generic. 

The set of all (weakly) generic mxm matrices with a given sign pattern is open and dense 
in the set of all mxm matrices with that sign pattern. 

Lemma 3.11. The rank r of a generic matrix A with connected graph G{A) is equal to the 
minimum of the number of rows or of columns of A. If G{A) has £ components Gi, . . . , Ge 
and is the minimal number of column or row nodes in Gi, then r = Yli=i ''hi- 



proof. Obvious. ■ 

Definition 3.12. For S G M''^'^' let S^ed denote a reduced iS-matrix, i.e., a matrix obtained 
from S by removing one column out of every pair of columns which are reverses of each other. 
Clearly, Sj-^d contains no reversible columns. The reduced flux pattern Uj-cd is obtained 
from S and S^cd- it is built from the sign pattern of —Sj^^ by setting all entries coming 
from positive entries in columns nonreversible in to 0. In particular, if all columns of S 
are reversible, then [/red is the sign pattern of —S^^,^. If no column of S is reversible, then 
Srcd = S and U^^d = U. 



Example 3.13. Let us revisit Example 13.51 A reduced S'-matrix S^cd and the reduced flux 
pattern U^cd are 



red 



an 





-ai3 


0,21 










^32 


033 


-a4i 











-^52 






red 



-[/ll 


-U12 





[^14 








-U22 


-[/23 





U25 


^31 















In most cases 5'rcd will be generic and hence of rank 3. By a straightforward computation. 



'S'rcd Uvcd 



-aiiUii - aisUsi 
-021 [/n 
OssUsi 





-auUi2 
-a2iUi2 — CI22U22 

^^32^^22 
041^^21 
(^52^^22 





— a22^723 

— a32^723 



0^52 1^23 



a2iUi4 


— a4i[/i4 






^22^25 
^^32 ^25 


— 0'52U25 
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After a possible renaming of the free variables in U^cd, SU = S^cdUrcd- This is the key 
observationwe use in the next sections in order to count or estimate the number of anomalous 
signs in cd(S'). ■ 

Lemma 3.14. If S is a real d x d' matrix, if S'rcd is any reduced S -matrix and U^cd the 
corresponding reduced flux pattern, then 

SU = S'rcd f^red 

after a possible renaming of the free variables in [/red- 



Proof. Suppose first that S = [ S^cd —Srcd ] • The corresponding matrix U is of the form 
(3.10) U-- 



Uo 



The sign pattern of Uq is the same as that of —Sj-d- Furthermore, nonzero entries of Uq 
coincide with negative entries of S'rcd- Similarly, nonzero entries of Ui coincide with positive 
entries of S'rcd - 

Clearly, SU = S^^diUo — Ui) = S'rcd^red (after a possible renaming of the free variables in 

Urcd)- 

Let us now look at the general case, where some of the columns do not have reverses in 
S. We 'expand' S to S = [ S'rcd —S'rcd ] by adding reverses of nonreversible columns. We 
insert rows of zeros at the appropriate places in U. Again, we write 



U 



Uo 
Ui 



As before, nonzero entries of Uq correspond to negative entries of S'rcd- (Nonzero entries of 



Ui correspond to a subset of the set of all positive entries of S'rcd-) As Uq 
concludes the proof. ■ 



Ui = Urcd, this 



3.5. Counting anomalous signs when S'rcd is square. Now we give our main theorem 
for square reduced S'-matrices. The result is strong and effectively reduces the problem to 
the matrix and graph-theoretic test of §2.4[ 

Theorem 3.15. Let S be a real dx d' matrix and suppose S'rcd "is a generic square invertible 
matrix. Then: 

(1) The number of terms in the core determinant cd(S') equals the number of terms in 
the determinant expansion o/t/red- 

(2) The number of anomalous signs of the core determinant cd(S') is the number of 
anomalous signs in the determinant expansion of Urcd- 

Remark 3.16. Note that the theorem gives a count of positive and negative terms in cd(S') 
when combined with Theorem 12.91 The number of (anomalous) signs in the determinant 
expansion of [/red is bounded above by the number of (anomalous) signs in the determinant 
expansion of the sign pattern of S'rcd- • 



Proof of Theorem 13.151 Follows immediately from Lemma 13.141 ■ 
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3.6. Rectangular S'rcd matrices. This section gives results and examples for the case of 
rectangular reduced S'-matrices. Our theorem for complicated situations would not easily 
yield the precise count. On the other hand, it yields estimates and in various simple cases it 
is effective. 

We can use Lemma fS.MI to provide a Binet-Cauchy expansion with fewer terms than there 
were in Lemma [3.71 namely: 

Lemma 3.17. For S G M'^^'^' having rank r, the core determinant is given by 
(3.11) cd(5) = (-1)'^-^ Yl det(^,ed(a|/?)) det{UUP\o^))- 

\a\,\f3\=r 

Theorem 13.151 and Remark 13.161 tell us how to count the number of positive, negative or 
anomalous signs in cd(S') with generic S'red- By the Binet-Cauchy formula ( 13. lip given in 
Lemma r3.17l we count the number of positive and negative terms for each of the det{Uj-cd{P\<y)) 
and take into account the sign of det{Sj-Ed{oi\P)). The sum of these will give us a count for 
the number of positive and negative terms in cd(S'). Note: due to the freeness of entries of 
f/, there is no cancellation between the summands. In particular, this count gives us a lower 
bound and upper bound on the number of anomalous signs in cd(iS'). 

Theorem 3.18. Suppose S E W^^'^' has rank r. Let S'red be a reduced S-matrix and U-^cd the 
reduced flux pattern. Suppose that S^^d is generic. 

(1) The number of anomalous signs in cd(S') is at least 

J2 m(t/red(/9|«)) 

\a\,\f3\=r 

and at most 

J2 t(f/red(/5|a))-m(^rcd(/9|a)). 

\a\,\(3\=r 

(2) The number of terms of sign (—1)'^"-'^ in cd(S') is at least 

m(f/red(/9|«)) 

and at most 

J2 t(t/red(/?|a))-m([/red(/3|«)), 

5,ed("l/3)eAr 

where M is the set of all r x r submatrices S'red that are not SD. 

Proof. (1) follows from the explanation given above, so we consider (2). For a SD matrix 
So = S'red(a|/?), all terms in the determinant expansion of the sign pattern of Sq have the same 
sign. Hence the same holds true for Uq = f/red(/9|«) which is the sign pattern of —Sq with 
possibly some entries set to 0. If \a\ = \P\ = r, then the sign of det([/o) is or (—1)^ times 
the sign of det(S'o). Hence by Lemma [3. 171 a term of sign (—1)'^"^ in cd(S') cannot come from 
a r X r SD submatrix of S'red- To conclude the proof, note that given Si = S'red £ A/", the 
term det (S'red («|/^)) det(f/red(/5|tt)) will contribute at least min{m_([/j), m+(?7j)} = m{Ui) 
terms of sign (—1)'^"-'^ in cd(S') and at most max{m_(f/j), m^{Ui)} = t{Ui) — m{Ui) terms of 
sign (-1)'^-^ (Here := ^red(/5|a).) ■ 
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3.6.1. Examples. 

Example 3.19. Let S he a d x d' matrix of rank r with generic S^ed- 

Suppose Sj-ed has no e-cycle interlacing with respect to a perfect matching, then cd(5') has no 
anomalous signs. 

This we now demonstrate. By assumption and Theorem 12.91 any r x r submatrix 5*0 = 
'S'rcd(«|/9) of Sred IS SD. Then J\f = 0, so by Theorem 13.181 cd(S') will have no anomalous 
signs. 

Conversely, ifcd{S) has no anomalous signs, then ^(t/rcd) has no e-cycles interlacing with 
respect to a perfect matching. 

To see why this is true, we invoke Theorem 12. 9[ Such an e-cycle and the perfect matching in 
G{Urcd) pick out a r X r submatrix Urcd{f3\(y) of Ured- The corresponding summand in 13.111 
will then yield at least one anomalous sign by Theorem | 



In the fully reversible case this yields a necessary and sufficient condition for cd{S) to have 
no anomalous signs. ■ 

Remark 3.20. We recall that for cfd(S') what we have just done is known [CF06] in the fully 
reversible case S = [ S^cd —S^cd ]■ What is shown in [CF06], implies that cfd(S') has no 
anomalous signs iff G(>5'red) has no e-cycles. ■ 

Example 3.21. Suppose 5'rcd is generic, S = [ S^cd —Srcd ] and the graph G{Srcd) is: 



I 



C2 



CI 



@ 



C3 



Cn 



There are n rows and n columns, so the rank of S (and Sj-cd) is n. G(S'red) supports exactly 
one rank n square matrix, Sj-^d itself. 

G{Sred) has one cycle with no c-pairs, so it is an e-cycle. G{Sred) admits 2 perfect match- 
ings: the e-cycle interlaces both matchings. So by Theorem 12. 9 j we get M is Sj-^d- 

Theorem 13.181 together with Theorem 12.91 imply 

1 = m(/7red) = Yl "^(^^) ^ "^('^) ^ Yl [^(^) - "^(^^)] = ^(^-d) - m([/,ed) =2-1 = 1. 

Thus generically cd(iS') has one anomalous sign (independent of t?, > 2). 

Alternative to Theorem 13. 18[ since S'red is square, we could have used Theorems 13.151 and 
12.91 which tell us that cd(S') will have 2 terms, one with a positive and one with a negative 
sign. 
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On the other hand, the number of anomalous signs in cfd(S') increases rapidly with n. 





^ of anomalous 


n 


signs in cfd(S') 


2 


1 


3 


2 


4 


5 


5 


13 


6 


34 


7 


89 


8 


233 


9 


610 


10 


1597 



This data is consistent with 

number of anomalous signs = Fib(2n — 3) 



(see the website http: //www. research. att.com/^njas/sequences/). We leave it to the 
interested reader to see if this is true. ■ 



Example 3.22. Suppose S 



[ Srcd -Srcd ] , -^red is gcueric and G(5'red) is the graph: 
R8) 



C4 



(R2) 
I 

I 



C2 



CI 



(g) 
I 



C3 



There are 8 rows and 4 columns, so rank of S and S'rcd is 4. One e-cycle. 

G'(S'rcd) admits 3 ■ 3 ■ 2 = 18 perfect matchings and the e-cycle is interlacing with respect 
to every one of those. Each perfect matching selects a 4 x 4 submatrix of S'rcd (or t/rcd) with 
one e-cycle in its graph. In total there are 9 such submatrices (each in M) with the graph 
of each one admitting two perfect matchings. 



Theorems 13.181 plus 12.91 imply 

9 = < m{S) < J2 [tiUi. 

Thus generically cd(5') has 9 anomalous signs. ■ 



m{U,)] = 9(2 



3.7. Few Anomalous Signs - An Algorithm. We have just looked at bounds for the 
number of anomalous signs in cd(S') for generic S-^ca- A small number of anomalous signs 
in the core determinant can be handled precisely using an algorithm we now describe which 
obtains necessary and sufficient conditions for cd(S') to have {zero or) one anomalous sign. 
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3.7.1. The zero-one anomalous sign algorithm: Suppose S is a. d x d' matrix of rank 
r. In order for the algorithm to work with certainty, we assume S^cd is generic. Let A/" be 
the set of all r x r submatrices of Sj-^d that are not SD. Given Si G J\f we use Ui to denote 
the corresponding submatrix of Uj-^d- We present the algorithm only for the case when cd(5') 
has no anomalous signs or the anomalous sign is 

Case E: Af has elements. 

Then cd(5') has no anomalous signs. 
Case N: A/" is nonempty. 

Subcase (a): All the Ui corresponding to Si & are SD. 

Take det(S'j) det{Ui) and look at its sign. If for all Si & Af this sign is (—1)'', then 
cd(S') has no anomalous signs. Otherwise for some Si & Af the sign is (— 
and the corresponding term det(S'j) det([/j) contributes t{Ui) terms with sign 
(—1)^"^ to cd(5'). If t{Ui) > 1, then there is more than one anomalous sign in 
cd{S). If there is Sj 7^ Si with sign(det(S'j) det{Uj)) = {-lY-\ then cd{S) will 
have more than one anomalous sign. Otherwise cd(S') has one anomalous sign. 

Subcase (b): There is exactly one Sq E J\f for which the corresponding Uq is not SD. 

If there is Si E Af \ {So} with the sign of det(S'i) det{Ui) equal to (—1)^, then 
cd(S') will have more than one anomalous sign. Otherwise we use the det sign 
test (Theorem 12.91) to compute m{Uo). 

(i) If m{Uo) > 1, then cd(S') will have more than one anomalous sign. 

(ii) Suppose m{Uo) = 1. If the number of terms t in det{Uo) is two, cd(S') will 
have one anomalous sign. So suppose t > 2. Let 

f +1 I m{Uo)=m+{Uo) 
1^—1 I otherwise. 

Now cd(S') will have one anomalous sign iff 

(3.12) esigndet(5o) = (-1)'^"^ 

(If f l3.12p fails, cd(5') will have more than one anomalous sign.) 
Subcase (c): There are at least two Si E Af for which the corresponding Ui is not SD. 
In this case cd(S') will have at least two anomalous signs. 

Lemma 3.23. The zero-one anomalous sign algorithm computes whether or not there is one 
{respectively no) anomalous sign. 

Proof. Case E is given in Example 13.191 Case N.(a) follows directly from the Binet-Cauchy 
formula (13.111) . For Case N.(b).(i), det(S'o) det([/o) has more than one anomalous sign, so 
cd(S') will have more than one anomalous sign. The proof of Case N.(b).(ii) is essentially 
contained in the statement. Finally, in the Case N.(c) two different Si contribute two different 
terms to the Binet-Cauchy expansion 13.111 for cd(S') each having at least one anomalous sign. 
■ 

Remark 3.24. The algorithm simplifies considerably in the fully reversible CcLSG, clS then Ui 
is SD iff Si is. Thus Case N.(a) cannot arise. Subcase (b) is equivalent to Af having exactly 
one element and Subcase (c) is equivalent to Af containing at least two elements. ■ 



This assumption is made purely for convenience of exposition. In fact, if Sred has at least two SNS r x r 
submatrices with nonsingular corresponding submatrices in C/rod, then this will automatically be the case. 



22 



J. WILLIAM HELTON, IGOR KLEP, AND RAUL GOMEZ 



Example 3.25. Let S = [ Sj-cd —Sred ], where 

' a 1 



1 

has rank 3 and a G M>o. Suppose a 7^ 2; this makes Sj-cd generic. The graph G{Sj.cd) is given 
by the following: 




The cycle Rl — CI— R2— C2— Rl has two c-pairs and is an e-cycle; it is the only e-cycle and 
it interlaces two perfect matchings. Both leave out R4 and select the same 3x3 submatrix 
5*0 of S'rcd- Hence Af = {So}. 

To count the number of anomalous signs in cd(S') apply the Algorithm l3.7.1l Our situation 
corresponds to Case N.(b) and we compute m(f/o), where Uq is the 3x3 submatrix of f/red 
corresponding to 5*0. By the det sign test, m{Uo) = 1 and t{Uo) = 3. It is easy to see 
that m{Uo) = m^lUo); thus by the zero-one anomalous sign algorithm, cd(S') will have one 
anomalous sign iff a — 2 = det(S'o) < 0. (If a > 2, cd(S') will have two anomalous signs.) ■ 



Another class of examples is cycles with short hair, a notion we now elucidate. A subgraph 
r of G is said to have short hair provided when edges of F are removed from G all paths 
in the remaining graph starting from a vertex in F have length < 1. 

Proposition 3.26. Suppose S'rod is generic and the graph G(5'rcd) is connected and contains 
at most one cycle and possibly some short hair {e.g. Example 13.51 or Example [CF05, Table 
1.1. (iii)]). Then the number of anomalous signs in cd(S') is < 1. 



Proof. Without loss of generality, G'(<S'red) contains a cycle S. If S is not an e-cycle, then we 
are in Case E of the algorithm and there are no anomalous signs in cd(5'). Thus we assume 
that S is an e-cycle. 

By Lemma [3.111 the rank r of Srcd is the minimal number of rows or of columns in S'^d- 
If r is bigger than the number of columns appearing in the cycle, then we are in Case E of 
the algorithm because any perfect matching will include some edge not in the cycle, thus 
making S not interlace it. Hence cd(S') has no anomalous signs. 

Otherwise r equals the number of columns appearing in the cycle. Then TV has only 
one element 5*0. The corresponding submatrix Uq of ?7red is either SD or its determinant 
expansion has two terms of opposite sign. Now the result follows from Case N.(a). ■ 



Note that results on the core determinant cd(5') given in this section have parallels for 
the Craciun-Feinberg determinant expansion cfd(S') which are easy to work out using the 
techniques in our paper. 
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4. Chemical Motivation 

This matrix theory paper is not directly aimed at producing chemical results but was 
inspired as an extension of the striking work of Craciun and Feinberg. We hope these 
extensions might someday prove valuable on chemical network problems and some methods 
they combine with are described in [CHWprept] and a consequence is Theorem 14.11 below. 

Now we turn to describing the connection between the core determinant from ^ and 
chemistry. 

A chemical reactor can be thought of as a tank with each chemical species flowing in 
(assume at a constant rate) and each species flowing out (assume in proportion to its con- 
centration in the tank). If the reaction inside the tank satisfles ^ = g{x), then when there 
are inflows and outflows, the total reaction satisfles 

$ = f{x) = g{x) + exin - 5x. 
at 

The Craciun-Feinberg determinant is the determinant of the Jacobian /' when 5 is 1 and 
it bears on counting the number of equilibria for this differential equation, cf. [CF05,CF06, 
CHWprept]. There is some discussion of small outflows vs. no outflows in [CF06iee]. 

The core determinant bears on a different problem. Assume the differential equation has 
reaction form f{x) = Sv{x). Let R (respectively i?^) denote the range of S (respectively 
its orthogonal complement); R is typically called the stoichiometric subspace. Let P be the 
projection onto R and P-*" onto R^. With no inflows and outflows, P-^f{x) = and clearly 
this implies the solution x{t) to the differential equation propagates on the affine subspace 

(4.1) := {x I P^x{t) = const = P^x°}. 

This reflects quantities (like the number of carbon atoms) being conserved. The flow on 

AixO has dynamics = '^*''^^"^f ^ ^ = Pf{Px + P^x^). Proposition 13.31 implies that the 
determinant of the Jacobian of this dynamics is the core determinant which we studied in 
this paper, namely, for any ^ in A^^-o 

(4.2) cd(5)(0 = det(P/'(0P). 

When cd(S') has no anomalous signs the degree theory arguments in §3 of [CHWprept] 
give a strong result for numbers of equilibria of the differential equation. 

Theorem 4.1. Suppose ^ = fb{x) := Sv^{x) has reaction form with v^{x) once continuously 
differentiable in x and depending continuously on a parameter < b < 1. Suppose each 
component Vj{x) of v''{x) is monotone nondecreasing. Suppose TW^o is compact. Suppose 
cd(<S') has no anomalous signs. 

If there are no zeroes fb{x) = for any b and any x on the boundary of Ai^o, then the 
number of zeroes for f^ in the interior o/A^^-o is independent ofb. 

The hypothesis that cd(5') has no anomalous signs can be weakened to cd(5')(^) does not 
equal for any ^ in A^^-o. 
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